Can randomness alone tune the fractal dimension? 
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We present a generalized stochastic Cantor set by means of a simple cut and delete process and 
discuss the self-similar properties of the arising geometric structure. To increase the flexibility of 
the model, two free parameters, m and b, are introduced which tune the relative strength of the two 
processes and the degree of randomness respectively. In doing so, we have identified a new set with 
a wide spectrum of subsets produced by tuning either m or b. Measuring the size of the resulting 
set in terms of fractal dimension, we show that the fractal dimension increases with increasing order 
and reaches its maximum value when the randomness is completely ceased. 
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I. INTRODUCTION 

The history of describing natural objects by geometry 
is as old as science itself. Traditionally, this has involved 
Euclidean geometry that restricts the natural bound- 
ary of our imagination to lines, circles, cones, spheres, 
cuboids and so on. However nature is not restricted 
to Euclidean shapes only. Instead, most of the natu- 
ral objects we see around us are so irregular and com- 
plex in shape that they can be described as geometrically 
chaotic, since they are not just merely complex but often 
contain different degrees of complexity. In 1975, Mandel- 
brot introduced the concept of fractal geometry to char- 
acterize these geometric monsters quantitatively. This 
extended our horizon of understanding and enables to ap- 
preciate that there exists some kind of order even in these 
seemingly complex and disordered geometric structural 
patterns. Often fractals are perceived as some brightly- 
colored computer generated patterns. This is due to the 
lack of enough analytically solvable theoretical models. 
The coastline of an island, a river network, the structure 
of a cabbage or broccoli, or even the networks of nerves 
and blood vessels in the normal human retina can be best 
described as fractals. Yet, more than twenty five years 
after they were first introduced, there is no generally- 
accepted definition of a fractal. However, it is loosely 
defined as a shape made of parts similar to the whole, 
in some sense. It is typically quantified by a non-integer 
exponent called the fractal dimension that can uniquely 
characterize the geometric structure. This definition im- 
mediately confirms the existence of scale invariance, that 
is, the objects look the same on different scales of obser- 
vation. To understand fractals, their physical origin and 
how they appear in nature we need to be able to model 
them theoretically. This is the main motivation of our 
present work. 

The simplest way to construct an ideal fractal is by 
repeating a given operation over and over again deter- 
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ministically. The classical Cantor set is such a simple 
text book example which is typically used for learning 
about the self-similar properties of an ideal fractal. It is 
created by dividing a line into n equal pieces, removing 
(n — m) of the parts created and repeating the process 
with the m remaining pieces ad infinitum |l|]. However, 
it is well understood that in our world nothing is sta- 
tionary or strictly deterministic and all natural complex 
objects, even intelligence and life, occur through some 
kind of evolution where randomness is an essential in- 
gredient. That is, nature favors randomness and fractals 
in nature appear through a continuous kinetic process, 
whereas the classical Cantor set is discrete in time. The 
simplest way to incorporate these essential ingredients is 
by generalizing the Cantor set. For example, instead of 
dividing an interval into n equal pieces, one can divide it 
randomly into n pieces and throw away (n — ni) of these 
also randomly 

The notion of random fractal has been widely used 
in almost every disciplines of science both theoretically 
and experimentally since its inception. Yet, the mecha- 
nism by which nature creates fractals and the relation- 
ship between the degree of order and the fractal dimen- 
sion is poorly understood. In this article, we introduce a 
stochastic process that can generate a continuous spec- 
trum of fractals (ranging from random to nonrandom) 
controlled by some intrinsic tuning parameters. This 
may be considered as a natural kinetic counterpart of the 
classical Cantor construction and which can be a poten- 
tial candidate in order to understand essential governing 
rules of creating complex geometric objects. We present 
a generalized version of the classical Cantor set. The 
rules of the process can be described as follows. The pro- 
cess starts with an initiator of the unit interval [0,1] and 
a generator which divides this interval into two pieces 
randomly, deleting some parts from either sides of both 
pieces at each time step. We introduce two intrinsic pa- 
rameters to determine the degree of randomness and the 
rate at which a given operation is to be repeated in order 
to create a fractal. 

The construction of this generalized Cantor set is not 
at all pedagogical. One immediate and potential applica- 
tion of the present model is the kinetics of the irreversible 
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and sequential breakup of particles that occurs in a vari- 
ety of physical processes and which has important appli- 
cations in science and technology. These include erosion 
grinding and crushing of solids [Q, polymer degra- 
dation and fiber length reduction j^, breakup of liquid 
droplets etc. to name just a few. In recent years, 
there has been an increasing interest in studying fragmen- 
tation, allowing variations that increase the flexibility 
of the theory in matching conditions of real phenomena 
such as the extension to higher dimension [Q, agglom- 
erate erosion mass loss fl , volume change ||T^ and 
fragmentation-annihilation UM- The generalized Cantor 
construction, as we formulated it, can describe the kinet- 
ics of fragmentation with a continuous mass loss [|[ ^. 
This is relevant in the kinetics of fragmentation of par- 
ticles where mass loss might occur due to evaporation, 
oxidation, sublimation, dissolution, melting etc., or even 
in the Yule-Furry process of cosmic shower theory where 
energy loss occurs due to collision Q. 

Another potential application of the present work is 
the first order phase transformation via the kinetics of 
nucleation and growth processes. The random deposi- 
tion of cuts in the generalized Cantor set can be thought 
of as the random nucleation of size-less seeds on the sub- 
strate (metastable phase). On the other hand, the con- 
tinuous removal of parts is in fact equivalent to the de- 
crease of substrate due to the growth of seeds (stable 
phase). The nucleation of a stable phase usually occurs 
due to large thermal fluctuations in the metastable phase 
and the subsequent growth occurs due to the competi- 
tion between the surface tension and the free energy den- 
sity difference between the stable and metastable phase. 
The classical work of such processes has been done in- 
dependently by Kolmogorov, Johnson, Mehl and Avrami 
(KJMA model) and has been extensively used to de- 
scribe non-equilibrium phenomena during phase trans- 
formation . In this model the nucleation of a stable 
phase, within the metastable phase, occurs at random 
obeying the Poisson stochastic process in space. In other 
words nucleation centers are homogeneously distributed 
throughout the entire metastable phase. Once nucleated, 
the stable phase keeps growing isotropically until it comes 
into contact with another similarly growing stable phase 
thus forming an interface. 

Interesting questions arise from the present work: (i) 
What is the role of the fractal dimension during the pat- 
tern formation? (ii) Is there any relation between the 
fractal dimension and the degree of order? (iii) What 
are the relevant parameters to tune the degree of order 
and what are their physical meaning? The present work 
is an attempt to answer these questions. 



II. FORMULATION OF THE CUT AND 
DELETE PROCESS 

In order to create a complex geometric objects here 
we give the general scheme for the formulation of the 



cut and delete model. To do so, we define the interval 
size distribution function or the concentration c(x, t) so 
that c(x, t)dx describes the number of intervals of size 
X at time t in the interval range x and x -t- dx. The 
distribution function c{x, t) then must have the following 
form of evolution equation 



dc{x,t) dc{x,t) 



dt 



dt 



dc{x, t) 



dt 



(1) 



delete 



The first term on the right hand side represents the ir- 
reversible and sequential cut process which can be de- 
scribed by the following master equation 



dc[x, t) 



dt 



= -c(x,t) / dyF{y,x-y) 

cut Jo 



dyciy,t)F{x,y-x). (2) 



Here F{x,y) is the breakup kernel describing the rules 
and the rate {a{x) = dyF{y,x — y)) at which an in- 
terval of size X + y breaks into sizes x and y. The two 
terms in Eq. (2) describe the destruction and creation of 
interval of size x due to the deposition of cuts on size x 
and larger than x respectively. The second term of Eq. 
(1) on the other hand, represents the deletion of parts 
of an interval from either ends soon after a cut is de- 
posited. The deletion of parts continues until these parts 
encounter another cut or the interval becomes a dust like 
interval; thereby stopping the loss of their masses. This 
deletion process can be expressed in the following differ- 
ential form 



9c(x, t) 



dt 



d{q{x)c{x,t)) 



delete 



dx 



(3) 



where q{x) describes the rate at which an interval of size 
X is lost by the deletion process. 



A. Model I: Random scission model 

We first consider here the simplest possible model of 
this kind, namely 



F{x,y) = 1. 



(4) 



This is typically known as the random scission model 
as it assumes that every points of the interval are 
equally likely to break. That is, the present model de- 
scribes the following recursive cut and delete process. We 
start with an interval of size x + y. At each time step, 
we choose a cut point with a uniform probability den- 
sity and produce two smaller intervals of size x and y. 
As soon as the new intervals are born, they start shrink- 
ing from either side, although the mechanism whereby 
they shrink would depend on the choice of q{x). In gen- 
eral, q(x) is assumed to follow a power law form with 
respect to particle size i.e. q{x) oc x'^ . However, once 
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the choice for the breakup kernel is made, the exponent 
7 is strictly restrictive, due to dimensional consistency. 
Therefore, in the present context it is essential to choose 
q{x) — mx^ where m is a positive, real and dimcnsionless 
constant. The rate equation for the distribution function 
c(x, t) then becomes 



dc{x, t) 
dt 



xc{x,t)-\-2 I c{y,t)dy + m 



dx'^c{x, t) 
dx 



(5) 



Instead of trying to solve the rate equation directly, we 
attempt to solve the rate equation for M„(t), the rfi^ 
moment of c(x,t), defined as M„(t) = x"c{x,t)dx 
with n > 0. Incorporating this definition into Eq. (5), 
we obtain 



dt 



+ n{l + 1/m) — 1/m 



f Mn+l{t). (6) 

(n + 1) 

Note that for m — the total mass or the first moment 
Mi{t) is a conserved quantity. However, for m > the 
system violates the mass conservation due to the contin- 
uous mass removal by the deletion process. The inter- 
esting feature of the above equation is that for m > 0, 
there are infinitely many conserved quantities depending 
on m. We can find these conserved quantities by simply 
inspecting the value of n = df for which the moment 
Mdf (t) is independent of time which is obtained by set- 
ting dMn{t)/dt = 0. The problem then rests on solving 
a quadratic equation in n. The real positive root of the 
equation then is 



(7) 



while the other root is a = —{df + 1 + m) a < 0. 
We can re- write the rate equation for Mn{t) as 



dMnjt) 

dt 



{n — df){n ^ a) 



A/„+i(t). 



(8) 



(n + l) 

Iterating the above equation enables us to obtain all 
the derivatives of Mn{t) at t = 0. Using this and the 
monodisperse initial condition into the Taylor series ex- 
pansion of Mn{t) about t — Q gives the explicit and exact 
solution for Mn{t) 



Mn{t) — 2Fi{n ~ df,n — a;n + I; —t) 



(9) 



where 2^1 (a, c; d) is the hypergeometric function. In 
principle, finding a proper integral transformation of Eq. 
(9) and comparing it with the definition of M„(t) would 
enable us to recover the exact solution for the distribution 
function c(x, t) (widely known as Charlesby's method 
[^6|). However, in the present work we are more in- 
terested in Mn(t) and its asymptotic temporal behavior 
than in c{x,t). The asymptotic expansion of the hyper- 
geometric function Q reveals that the exact time de- 
pendence of Mn{t) is Mn{t) The asymptotic 
decay of M„ (t) is thus linear in n which implies that the 
system obeys a simple scaling. This will be shown vigor- 
ously in the following section. 



1. Scaling theory 

A closer inspection of Eq. (5) immediately reveals that 
there are two governing parameters x and t that com- 
pletely describe the governed parameter c{x,t). This is 
due to the fact that the dimension of the two variables x 
and t are interlocked by the choice of the breakup kernel 
and through the rate equation itself. That is, x, the co- 
factor of c{x, t) of the first term on the right hand side of 
Eq. (5), must bear the dimension inverse of time t. This 
enables us to choose any of the two as an independent pa- 
rameter while the other can be expressed in terms of this. 
That is, we can define ^ = te as a dimcnsionless quantity. 
So, if we assume t to be the independent quantity then 
c(x, t) can be expressed in terms of t alone i.e. 



(10) 



Here 6 is yet an undetermined exponent and must take 
the value for which t^ bears the dimension of concen- 
tration or the distribution function c. This means that 
all plots of vs xt for any initial condition must fall 
onto each other which is the hallmark for the existence 
of scaling typically known as data collapse formalism [ p^ . 
Substituting the scaling ansatz (Eq. (10)) into the defi- 
nition of Mn{t) and demanding that Md^ be a conserved 
quantity, one immediately obtains 9 = 1 + df. 



2. Fractal analysis 




FIG. 1: Fractal dimension versus mass removal strength m. 

Associating broken objects with disorder is perhaps 
the most immediate and natural thought that cross our 
mind. Indeed, an assembly of fragmented intervals can 
be considered as one of the simplest yet most common 
example of a disordered system. The existence of scaling 
is already a proof that an assembly of such sequentially 
broken intervals has some degree of regularity. Owing 
to the random nature of the present model and due to 
the mass removal term, it is clear that when the process 
evolves ad infinitum, it will not create the whole set that 
describes a line. Rather, it will produce a subset of the 
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full set (line) . In order to measure the size of the set cre- 
ated in the long time limit, we define the average interval 



^^^^ ^ ir ^'^(^' ^)^^ _ ^-1 
c(x, t)dx 



(11) 



This shows that the mean or typical interval size de- 
creases in time algebraically. We can use this as a yard 
stick and count the number of the segments we need to 
cover the resulting set. We find that the number of yard 
stick N{6) scales as 



lim N(S) 
s — >a 



(12) 



where N = c{x,t)dx. The exponent df is known 
as the fractal dimension or the Hausdorff-Besicovitch di- 
mension of the resulting set created by the cut and delete 
model. The distance between the points in space is the 
key to the definition of the Hausdorff-Besicovitch dimen- 
sion ijl^. Note that d/ is the root of the quadratic equa- 
tion in n obtained by setting '^^^J^ '^^ = (a condition 
for stationary solution). The expression for dy, Eq. (7), 
reveals that as the value of m increases, the fractal dimen- 
sion decreases very sharply and in the limit m — > oo, 
df — > 0. Fig. (1) does indeed show that as m increases 
the strength of the mass deletion term increases with 
respect to the terms responsible for describing the cut 
processes. This means that as m increases, the size of 
the arising set decreases sharply due to the fast disap- 
pearance of its member. While as m — > 0, d/ — > 1, 
that is, we recover the full set from which the subset is 
derived. 



B. Model II: Gaussian cut and delete model 

Next we treat the following quadratic type of breakup 
kernel 



F{x,y) = (xy) 



b-l 



(13) 



This is a generalized version of the previous model and it 
describes the following recursive cut and delete process. 
We start with an interval (x -I- y) and at each time event 
we choose a cut point with a Gaussian probability density 
if & > 1. The spread of the density of the cut around 
the center of the interval depends on the values of b as 
we will show later. Of course we recover the previous 
random scission model when b — 1. Our principal goal of 
studying this model is to understand the role of b during 
the cut and delete process. In order to extract the role 
of 6, it is essential to check whether the relative strength 
between the cut and delete process changes as we change 
the b value. Note that the breakup rate a{x) for this 
model is a{x) = px'^''-^ where p = [T{b)f /T{2b). It is 
thus clear that as b increases so does the relative strength 
of the cut process with respect to the delete process if we 
do not tune the m value accordingly. That is, to judge 



the role of b we ought to give equal weight to all the terms 
of Eq. (1) so that each of them can compete on an equal 
footing. This can only be done if we set q{x) — px^'' 
(i.e. m = p) and hence the relative strength between the 
two terms of Eq. (1) stays the same as the value of b 
increases. The rate equation for Mn{t) then is 



dMn 
dt 



= -m 



r{2b) 



(n+l)-^.^jAW_,. 

(14) 



The problem now rests on finding the n value for which 
the moment becomes a time independent quantity. From 
the previous fractal analysis we learnt that the real posi- 
tive root of the polynomial equation (obtained by setting 
dMn (t) /dt = {]) is in fact the fractal dimension of the re- 
sulting set. Applying the scaling theory based on dimen- 
sional analysis, like we did in the previous section, gives 
9= (l+(i/)z where the kinetic exponent is z = 1/(26—1). 
A detailed numerical survey reveals that the fractal di- 
mension increases monotonically with increasing b (see 
Fig. (2)). In order to find the fractal dimension in the 
limit b — > oo, we can use the Stir lings approximation in 
the resulting polynomial equation. In doing so we obtain 
the following equation 



ln[n + 1] + (n - 1) ln[2] = 



(15) 



and find df — 0.45699956 is the only real positive root 
and hence is the fractal dimension df. In an attempt to 
search for the role of b we study the following model 



F{x, y) ^ {x + y)'^S{x - y). 



(16) 



It describes that at each time an interval is chosen and 
the cut is placed exactly in the middle to produce two 
equal sized smaller intervals. It is note worthy to men- 
tion that the homogeneity index gamma only tunes the 
rate a{x) = x'^ /2 and does not play any role in control- 
ling the location of each cut on the interval vis-a-vis the 
fractal dimension. The exact value of 7 is not important 
provided 7 > and therefore we can choose 7 = 26 — 1 
without loss of generality. Once again the dimensional 
consistency requires that q{x) — x'^^/2 while the factor 
1 /2 is to account for the equal weight in the cut and delete 
process. Introducing this model into Eq. (1) and insert- 
ing the definition of moment we obtain the following rate 
equation for M„(t) 



dMnjt) 
dt 



1) 



-2b-l 



it). 



(17) 



In order to find the fractal dimension for this process, 
we set the numerical factor on the right-hand side of this 
equation equal to zero and then take the natural loga- 
rithm on both sides to obtain the n value for which M„ (t) 
is time independent. In doing so, we arrive at the same 
functional equation for n as we found for the Gaussian 
model with 6 — > 00, namely Eq. (15). 

This is in fact the main message of the present work. It 
shows that the kernel F{x,y) = (a;?/)^~^ behaves exactly 
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FIG. 2: Fractal dimension versus b. 



in the same fashion as for F{x, y) = {x + y)'^^~^5{x — y). 
We thus find that in the hmit b — > oo, the result- 
ing distribution of points is a set with fractal dimension 
Df = 0.45699956, which is a strictly self-similar fractal 
as randomness is ceased by dividing the intervals into 
equal pieces. We are now in a position to give a physical 
picture of the role played by b. First of all, the model 
with 6=1 describes that the frequency curve of placing 
cuts about the size of the fragmenting particles is Pois- 
sonian in nature. Consequently, the system enjoys the 
maximum randomness and the corresponding fractal di- 
mension is df = 0.414213. For 6 > 1, the frequency curve 
of placing cuts about the size of the fragmenting parti- 
cles is Gaussian in nature, meaning that as the value of 
b increases, the particles are increasingly more likely to 
break in the middle than on either end. That is, as b 
increases, the variance decreases in such a manner that 
in the limit b — > oo, the variance of the frequency curve 
becomes infinitely narrow, meaning a (5-f unction distri- 
bution for which the fragments are broken into two equal 
pieces. Therefore, there is a continuous spectrum of frac- 
tal dimensions between 6—^0 when df = 0.414213 and 
b — >oo when df = 0.45699956 (see Fig. (2)). 



III. DISCUSSION 

In this work we have presented a cut and delete process 
that can capture the notion of the Cantor set construc- 
tion. It has potential applications both in the kinetics 
of fragmentation and in the kinetics of nucleation and 
growth phenomena. We have found that in the long time 
limit the system creates a geometric pattern which is sta- 
tistically self-similar. To quantify the resulting geomet- 
ric structure, we invoke the idea of fractal geometry and 
show that the fractal dimension is strictly dependent on 
the strength of the mass removal term m. Thus by tun- 



ing this m, we obtain a fractal dimension of any value 
between < d/ < 1. We have further generalized our 
model allowing a parameter b that can tune the degree 
of randomness by controlling the way the cuts are de- 
posited. This is done in an attempt to find a precise 
answer of how fractal dimension changes as we tune the 
degree of order of the resulting geometric structure. We 
show that when 6=1 the system enjoys the maximum 
degree of freedom in the sense that cut points are cho- 
sen with uniform probability density and so is the delete 
process, since the two processes are interlocked via di- 
mensional consistency. For 6 > 1, at each time event a 
cut point is more likely to be chosen in the middle than 
on either end of a given interval. That is, the frequency 
curve of the cut density about any given interval size is 
Gaussian in nature and as 6 increases the middle point 
is increasingly more likely to be the cut point than any 
other points in the interval. The b value is intimately re- 
lated to the variance of the probability density of cut. As 
6 increases, the variance of the probability density curve 
decreases in such a manner that in the limit 6 — > oo, the 
variance of the curve reduces to an infinitesimally narrow 
line. The equivalence of the 5-function and the Gaus- 
sian model in the b — > oo limit does indeed prove this. 
Changing the 6 value effectively means we are creating 
an increasingly ordered set. In this work we have derived 
a spectrum of fractal dimensions between df ~ 0.414213 
for 6 = 1 and df = 0.45699956 when 6 — > oo. Fig. (2) 
and a detailed numerical survey confirms that the fractal 
dimension increases monotonically with increasing order 
and reaches its constant value when the randomness is 
completely ceased. 

In summary, we have discussed the scaling theory of 
a particular class of cut and delete process, emphasizing 
on the dimensional analysis. In order to make the scaling 
theory more meaningful than what a simple scaling can 
provide we have invoked the concept of fractal analysis. 
We have identified a new set with a wide range of sub- 
sets produced by tuning the degree of randomness only. 
We have quantified the size of the resulting set obtained 
in this way by fractal dimension and have shown that 
the fractal dimension increases with increasing order and 
reaches its maximum value when the pattern described 
by the set is perfectly ordered, which is contrary to some 
recently found results [Q. Finally, we argue on the ba- 
sis of our findings that fractal dimension, degree of order 
and the extent of ramifications of the arising geometric 
patterns are intimately connected to each other. 
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